Surface Critical Behavior in Systems with Absorbing States 
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We present a general scaling theory for the surface critical behavior of non-equilibrium systems 
with phase transitions into absorbing states. The theory allows for two independent surface expo- 
nents which satisfy generalized hyperscaling relations. As an application we study a generalized 
version of directed percolation with two absorbing states. We find two distinct surface universal- 
ity classes associated with inactive and reflective walls. Our results indicate that the exponents 
associated with these two surface universality classes are closely connected. 



PACS numbers: 05.40. +j, 64.60. Ht, 68.35.Rh 



The critical behavior of systems with boundaries has 
been the focus of much research in recent years 13] . So far 
most work on surface critical behavior and on the analysis 
of surface universality classes has been within the frame- 
work of equilibrium statistical mechanics. However, the 
same ideas and principles also apply to non-equilibrium 
systems. A prominent example of such a non-equilibrium 
process is directed percolation (DP), which is the generic 
model for systems with a non-equilibrium phase transi- 
tion from a state with activity (e.g., with a nonzero den- 
sity of particles) into a so-called absorbing state (with 
zero activity) . An understanding of DP is important for a 
wide variety of different systems encompassing epidemics, 
chemical reactions, interface pinning/depinning, spatio- 
temporal intermittency, the contact process, and certain 
cellular automata Recently, however, studies have 
been made of a number of systems with absorbing states 
which do not belong to the DP class. One prominent 
example is a particular reaction-diffusion model called 
branching and annihilating random walks with an even 
number of offspring (or BAW for short, where in this pa- 
er BAW refers exclusively to the even offspring case) 
Other systems in the BAW class (at least in 1 + 1 di- 
mensions) include certain probabilistic cellular automata 
monomer-dimcr models jsj, non-equilibrium kinetic 
Ising models H, and generalized DP with two absorbing 
states (DP2) |j. 

In this paper we address the impact of walls on sys- 
tems with phase transitions into absorbing states. We 
have developed a general scaling theory which allows for 
two independent surface exponents, which satisfy gener- 
alized hyperscaling relations. As an application, we have 
investigated the surface critical behavior of DP2. Our 
numerics indicate that DP2 exhibits a far richer surface 
structure than DP: we find two different surface univer- 
sality classes for DP2 with inactive and reflective walls, 
and our numerical results indicate that the exponents 
associated with these two classes are closely connected. 
These results can be successfully contained within our 
scaling theory. However, we emphasize that the theory 
is much more general than this and should also apply to 



other types of systems with walls and absorbing states, 
e.g., to surface effects in catalytic reactions and systems 
exhibiting self-organized criticality ||. 

Before turning to the surface critical behavior of DP2 
(in I + 1 dimensions) and BAW, we begin by discussing 
the main features of the corresponding bulk systems and 
then identify some differences and similarities with DP. 
Many models in the BAW class conserve particle 

number modulo 2, but this appears not to be the funda- 
mental requirement for the emergence of the new univer- 
sality class. Instead the key underlying feature seems to 
be the presence of a symmetry relating the various ab- 
sorbing states ||. This has been further demonstrated 
by Hinrichsen who introduced a generalized version of the 
Domany-Kinzel model with n absorbing states Q. This 
model, which we will refer to as DPn, is defined on a d- 
dimensional lattice (in space). At time t, the state s* of 
the i-th site can be either active (A) or in one of n inactive 
states (Ji, . . . , I n ). In 1 + 1 dimensions, the update prob- 



abilities Pfs' +1 |s* 



are given by P(I k \I k ,I k 



I, P(A\A, A) = l-nP(I k \A,A) = q, P{A\I k ,A) = 
P{A\A,I k ) = p, P(I k \I k ,A) = P(I k \A,I k ) = l-p, 
P(A\I k ,Ii) = I, where (k,l = 1, ...,n;k ^ /) (see also 
[Q for a more complete explanation of the model) . For 
n = I these rules are equivalent to the Domany-Kinzel 
model which belongs to the DP universality class (apart 
from one special point which belongs to the compact DP 
universality class) |l0|,[nj. For n > 2, the distinction 
between regions of different inactive states is preserved 
by demanding that they are separated by active ones. 
Monte Carlo simulations show that bulk DP2 belongs to 
the bulk BAW class in 1 + 1 dimensions |?J , whereas this 
probably does not hold in higher dimensions. 

The growth of both BAW and DP clusters in the bulk 
close to criticality can be summarized by a set of inde- 
pendent exponents. A natural choice is to consider v± 
and which describe the divergence of the correlation 
lengths in space, £j_ ~ |A| _l/i , and time £|| ~ |A| _i/ ii, 
where A = p — p c describes the deviation from critical- 
ity. We also need the order parameter exponent (3, which 
can be defined in two a priori different ways: it is either 
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governed by the percolation probability (the probability 
that a cluster grown from a finite seed never dies), 

P(A)~A A ~ d , A > 0, (1) 

or by the density of active sites in the steady state, 

n(A)~A**«», A>0. (2) 

For the case of DP, it is known that f3 is unique: /3 S eed = 
/3dons in any dimension. This follows from theoretical 
considerations |l2|]lg| l and has been verified by exten- 
sive numerical calculations. The relation also holds for 
BAW in 1 + 1 dimensions, a result first suggested by 
numerics and now backed up by an exact duality map- 
ping fl4| . However, this exponent equality is certainly 
not always true — for example it breaks down for certain 
systems with infinitely many absorbing states ||il||l6|| . 

Furthermore, /3 soe d ^ Pdcns for BAW in high enough 
dimension: if we consider the mean-field regime valid for 
spatial dimensions d > d c = 2, then the system is in 
an inactive state only for a zero branching rate, whereas 
any non-zero branching rate results in an active state. 
The steady-state density (||) approaches zero continu- 
ously (as the branching rate is reduced towards zero) 
with the mean-field exponent Pdens = 1- Nevertheless, 
for d > 2 the survival probability (0) of a particle cluster 
will be finite for any value of the branching rate, implying 
that Pseed = in mean-field theory. This result follows 
from the non- recurrence of random walks in d > 2. 

From the perspective of formulating field theories for 
BAW, the 1 + 1 dimensional case poses considerable diffi- 
culties [0 . These stem from the presence of two critical 
dimensions: d c = 2 (above which mean-field theory ap- 
plies) and d' c « 4/3 (where for d > d' c the branching 
reaction is a relevant process at the pure annihilation 
fixed point, whereas for d < d' c it is irrelevant there fin\). 
This means that the (physically interesting) spatial di- 
mension d = 1 cannot be accessed using controlled ex- 
pansions down from the upper critical dimension d c = 2. 
However if we assume that a (bulk) scaling theory can 
be properly justified (as it can be for DP, and BAW for 
d > d' c ), then it is straightforward to relate the above set 
of exponents to those of other quantities. Keeping the 
distinction between /3 see d and /3 d ens, the average lifetime 
of finite clusters, (t) ~ |A| _T , satisfies t = un — /3 sco d, 
and the average mass of finite clusters, 



(a) ~ |Ar , 
leads to the following hyperscaling relation: 
z/ii + dv± = (3 sccd + dens + 7. 



(3) 



(4) 



Note that (|J) is consistent with the distinct upper critical 
dimensions for BAW and DP. Using the above mean-field 
values for BAW and v\_ = 1/2, v\\ = 1, and 7 = 1, we 
verify d c = 2. In contrast, for DP one has the mean-field 
exponent /3 scc d = 1 and d c = 4. 



We now turn to the surface critical behavior of DP2 
and show how the above relations and exponents are 
modified in a semi-infinite geometry where we place a 
wall at x± = [x = (xy , x±), with the _L and || directions 
being relative to the wall]. In the simulations we start 
from an absorbing state, where all sites are in the state 
ii. We then initiate a cluster by placing a seed (site in 
state A) next to the wall. However, the analogy with DP 
is no longer immediate, as our numerical measurements 
in 1 + 1 dimensions indicate that DP2 supports an addi- 
tional surface exponent as well as an additional surface 
universality class. The type of surface universality class 
is governed by the choice of boundary condition (BC). 
We have studied two types of BC: the inactive BC (IBC) 
where the wall sites are always in the inactive state 7i, 
and the reflective BC (RBC), where the wall acts like a 
"mirror" by letting imaginary sites next to the outer side 
of the wall be the mirror images of those on the inside. 

By growing a DP cluster near an IBC wall, it has been 
observed numerically in d — 1,2 that certain exponents 
are altered [ ^8|Jl9|. This behavior has been explained by 
a scaling theorypOfl that explicitly takes surface critical 
phenomena into account and connects IBC with the ordi- 
nary transition pi] . Apart from the above (three) inde- 
pendent bulk exponents, an additional universal surface 
exponent must be included, which satisfies a generalized 
hyperscaling relation |20|. The survival probability (Q) 
for a cluster started close to the wall has the form 



(5) 



where the subscript '1' refers to the wall. However, in 
analogy with the bulk case, an order parameter can also 
be defined by the density of active sites on the wall in 
the steady state: 



"1 



(A) ~ A 01 ^ 



A > 0. 



(6) 



More generally the steady-state density (g) is now given 
by n(A,x±) = A^ dons <p(x±/£±), where the scaling func- 
tion ip behaves in such a way that n(A, x±) for x±/£±_ -C 
1 crosses over to the surface behavior (|J). 

For the case of DP, the surface exponents fulfill 
/3i,scod = /3i,dcns, as can be shown by a field-theoretic 
derivation of an appropriate correlation function [ p0| . 
However, for DP2 this exponent equality is no longer 
true. Our numerical results in 1 + 1 dimensions yield two 
distinct surface exponents, (3i, SC cd 7^ Pi, dens, although the 
corresponding bulk exponents coincide, as expected. The 
values of these surface exponents depend on the bound- 
ary conditions and by changing from IBC to RBC or vice 
versa, we observe that the assignment of the exponents 
is interchanged (see below). Further investigations are 
needed in order to determine whether the wall may have 
broken a (duality) symmetry present in the bulk (which 
forces the bulk exponents to coincide) and whether the 
operation of this symmetry relates IBC to RBC and vice 
versa. In contrast for surface DP, we note that IBC and 
RBC belong to the same surface universality class. 
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By keeping Pi. SC cd and /?i,dons distinct, we can now set 
up a general scaling theory for the surface critical behav- 
ior in systems with absorbing states. An ansatz for the 
coarse-grained density of active sites p\ at the point (x, 
t) of a cluster grown from a single seed located next to 
the wall, has the form 



Pi (x,t, A) = A*.— (x/U, 



(7) 



The A-prefactor comes from (||) for the probability that 
an infinite cluster can be grown from the seed, and from 
(||) for the (conditional) probability that the point (x, t) 
belongs to this cluster. The shape of the cluster is gov- 
erned by the scaling function f\ and we assume that the 
density is measured at a finite angle away from the wall. 
If the density is measured along the wall, we have instead 



Pn(x,t, A) = A*— /n (x/U, 



(8) 



as we pick up a factor Aft dens rather than A^ dons for 
the probability that (x, t) at the wall belongs to the 
infinite cluster. In 1 + 1 dimensions, (||) reduces to 
Pii(t > A) = A/ J i.—+ft.-»/ u (t/e|,). 

Starting from a seed on the wall, the average lifetime 
of finite clusters, (t) ~ |A| _Tl , satisfies n = ^ — /3i, se ed- 
The average size of finite clusters follows from integrating 
the cluster density (Q) over space and time: 



(9) 



where the surface (susceptibility) exponent 71 is related 
to the previously defined exponents via 



dv 1 



ft,. 



seed H~ Pdens ~t~ Tl ■ 



(10) 



The only difference from (||) is that we have now included 
a wall. Analogously, by integrating the cluster wall den- 
sity (||) over the (d — l)-dimensional wall and time, we 
obtain the average (finite) cluster size on the wall, 



(Swall) 



where 



(d-l) 



l.sccd 



l,do 



+ 7i,i- 



(11) 



(12) 



Note that if the 7 susceptibility exponents obtained from 
( |l0| ) and ( |l2] ) are negative, then they should be re- 
placed by zero in (^) and (]TT|) . The above scaling the- 
ory is generic since it allows for /3i, scc( j and /3i,dens to 
be independent surface exponents. When the scaling 
theory is applied to DP, it can be fully justified (with 
/3i, S ccd = /3i,dcns) @- However, if we apply the the- 
ory to BAW p^| , it would again be desirable to obtain 
a secure renormalization-group justification for the scal- 
ing behavior. In particular it would be important to de- 
termine from the field-theory whether two independent 
surface exponents are present. However given the funda- 
mental difficulties encountered already in the bulk field- 
theoretic analysis of BAW in 1 + 1 dimensions |L7| , this 



kind of analysis for the surface is unlikely to give a com- 
plete justification of the scaling theory. 

In order to confirm our scaling theory we have per- 
formed numerical simulations for DP2 in 1 + 1 dimen- 
sions with walls constrained by IBC or RBC (see |2^| for 
details). We have also performed simulations for DP2 
without a wall and obtained results for the exponents 
in complete agreement with |Q. There are several esti- 
mates for fldens (= /3 S ecd) available pi] : we have used 
/3 do „ s = 0.922(5) || . 

We extract the critical exponents from several mea- 
sured quantities. Using (J5J) , we find that the survival 
probability for the cluster to be alive at time t has the 
following behavior at criticality (A = 0) 



Pi(t) ~t 



#l,seed — Pl,seed/V\\ 



(13) 



Integrating the densities (^) and fll) gives expressions for 
the activity at criticality as function of time [p3[, e.g., 



Ar x (i)~r 



K l = d\ - ^dens - h 



seed; 



(14) 



where we have introduced the envelope (or "roughness" ) 
exponent \ = vx/vh, and <5dcns = Aiens/^il' Note that 
( |l4| ) corresponds to the hyperscaling relation ( |l0| ) at crit- 
icality with 71 = v\\{l + Ki), since (s) ~ J dtNi(t). For 
further confirmations of our numerical data we also con- 
sidered the cluster size distributions at criticality. The 
cluster size s scales as s ~ £^£iin(A) ~ A -1 / ', with 
1/(7 = dv± + z/|| — Pdens- The probability to have a clus- 
ter of size s then reads |23l 



Pi{ s ) 



-Ml 



Mi = 1 



Ml.seed 



dv±_ + V\\ - Pde 



(15) 



In Table [| we list our estimates for the critical expo- 
nents for DP2, where #i,dens = /^l.dcns/^n is obtained 
from (||) by measuring the activity at the wal l p3| , and 
M = 1+A;ccd/ {di'x+^W ^ /3dcns) corresponds to (|15|) in the 
absence of a wall. The results are in complete accordance 
with our theoretical analysis: bulk exponents are unal- 
tered whereas the wall introduces two separate surface 
exponents. We have also carried out bulk and surface 
simulations for A > and confirmed that our data could 
be collapsed according to an appropriate survival proba- 
bility scaling function [see @ for the surface case] , using 
our exponent estimates. This numerically confirms the 
validity of the relation 6 — for the bulk as well as 

for both sets of corresponding surface exponents [£7j . We 
further observe that the IBC, RBC boundary conditions 
lead to different exponents thus showing the existence 
of two distinct surface universality classes. Furthermore, 
/?i,seed 7^ /3i densj although by changing BCs we observe to 



good accuracy that f3y j^2 



R (RBC) R {RBC) _ R (IBC) 
"l.dens ' ^l jSe cd ^l.dens • 



As noted above, this suggests that the two BCs for DP2 
are related by a symmetry. By universality, we expect 
the same relations to apply to BAW |||| . 

By using the explicit definitions of IBC, RBC we can 
argue that /3i jS ccd and /3i dens should indeed depend on 
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the BCs. There will be more activity next to the wall for 
IBC than for RBC, since the latter can have regions of 1% 
located at the wall. Once created, these 1% regions will 
survive until the activity returns to the wall. Thus, from 
the wall density (§), it follows that < . On 

the other hand, the existence of these I2 regions implies 
that the survival probability for IBC will be smaller 
than for RBC, leading to p[ T ae ^ > Pi^d ■ However, 
from the observation that /3i iSee d + /3i,dens is indepen- 
dent of the BC, it follows that the average mass on the 
wall ([ll]) is the same for IBC and RBC. We have also 
studied several other BCs and found that these give the 
same scaling behavior as either RBC or IBC depending 
on whether or not the above-mentioned I2 regions can 
disappear only at the wall or also in the bulk pq] . In 
terms of the BAW model, however, the distinction be- 
tween the two BCs is slightly different: IBC respects the 
"parity" symmetry of the bulk, whereas RBC breaks it. 

For DP it has been customary to investigate whether 
the critical exponents can be fitted by (simple) ratio- 
nal numbers p6| . Such a fitting has also been tried for 
bulk BAW with the followingguesses in 1 + 1 dimensions: 
k — x — 2<5 = and \ = 4/7 || . These estimates lead im- 
mediately to S = 2/7 (and f3/v± = 1/2). It is intriguing 
to note that our numerical results for DP 2 in addition 
suggest that Hi =3/2 for IBC and 4/3 for RBC. From 
Eq. (|l5|), it then follows that 5i, sccd = 9/14 for IBC and 
3/7 for RBC. We would need one more relation in order 
to obtain the last independent exponent. In fact, we ob- 
serve numerically that 2v\\ — /3i !SCO d — /3i.dcns = 3 (28), is 
valid to within one percent [^9| . 

In conclusion, we have presented a generic scaling the- 
ory of surface critical behavior in systems with absorbing 
states. In particular we have for the first time studied 
the surface critical behavior of DP2, a model belonging 
to the BAW universality class in 1 + 1 dimensions. Nu- 
merical simulations of the DP2 model with two different 
types of boundary conditions have uncovered two sur- 
face universality classes. Our most important result is 
that two surface exponents are required to describe the 
surface critical behavior. The results also indicate that 
the exponents associated with these two surface univer- 
sality classes are closely connected. We emphasize that 
our theory is generic for systems with absorbing states 
and therefore should also apply to surface effects in, for 
example, systems exhibiting self-organized criticality. It 
would also be possible to generalize our theory to allow 
for edges and corners, which would introduce new expo- 
nents and other hyperscaling relations. 
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TABLE I. Critical exponents obtained from our simula- 
tions. For comparison we also list the exponents for DP with 
an IBC wall in the first column ^.jigg. 





DP (IBC) 


DP2 


DP2 (IBC) 


DP2 (RBC) 


^dens 
/3d ens 


0.159 47(3) 
0.276 49(4) 


0.287(5) 
0.922(5) 


0.288(2) 
0.93(1) 


0.291(4) 
0.94(2) 


^l,seed 
/?l,seed 


0.4235(3) 
0.7338(1) 




0.641(2) 
2.06(2) 


0.426(3) 
1.37(2) 


^l,dens 
jj\ ,dens 


0.4235(3) 
0.7338(1) 




0.415(3) 
1.34(2) 


0.635(2) 
2.04(2) 


M 
Mi 


1.108 25(2) 
1.2875(2) 


1.225(5) 


1.500(3) 


1.336(3) 
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